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Over the past few decades, following the discovery of chaotic behaviour in deterministic 
nonlinear dynamical systems, there has been considerable interest in the study of nonlinear 
differential equations and nonlinear difference equations with applications ranging across 
physics, chemistry, biology and economics. In many of these studies a knowledge of the 
properties of simple one-dimensional maps such as the logistic map, the dissipative standard 
map, the Lorentz map, the Henon map, the tent map, the quadratic map, etc., has provided 
fundamental insights[l — 5]. This has been particularly the case in understanding determin- 
istic chaos, understanding various transitions to chaos as a control parameter is varied, and 
identifying universality classes for the onset of chaos similar to universality in the theory of 
critical phenomena. 

Most detailed studies of one- dimensional maps have been confined to maps with a single 
control parameter and many of the chaotic and universal features of these maps are now 
well understood. For example, Feigenbaum[5] identified universal values for convergence 
ratios in the period-doubling route to chaos in certain one-dimensional maps. Universality 
has also been identified in the quasiperiodic route to chaos and in intermittency[3]. Two 
of the most useful quantitative measurements of dynamical properties of single parameter 
one-dimensional maps have been furnished by Lyapunov exponents and fractal dimensions. 
In this paper we present the results of our measurements of these and related quantities, 
multifractal spectrum, for a one-dimensional map with two control parameters - the Kim- 
Kong map [6]. The motivation for this work is twofold; firstly explore multifractal analysis 
for discrete dynamical systems and secondly to investigate the possibly competing effects of 
two control parameters in discrete dynamical systems. 

The majority of the applications of multifractal analysis[7 — 9] to date have been to 
physical systems. For example in recent applications, Kim and Kong used the box-counting 
method to compute generalized dimensions and scaling exponents for mountain heights[10] 
and sea-bottom depths [11] in Korea. 
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The Kim-Kong map is a two parameter map of the unit interval onto itself defined by [6] 

x(n + l) = /(ar(n), 7 ,/?) 

= 7 exp[-/3(logx(n)) 2 ][l-exp[- / 3(logx(n)) 2 ]] (1) 

with x{n) G (0, 1), where 7 and (3 are the two control parameters. To begin with we note 
that the Kim-Kong map is a dissipative unimodal map with a locally quadratic maximum. 
The maximum value is f(x) = ^ which occurs at x = exp( — y^J). A series expansion about 
the maximum yields 

/(*) « \ ~ 7/?exp(-2^?)(x - exp(-2y^)) 2 , (2) 

and the Schwarzian derivative is nee; ative for all x G (0, 1). If \ > exp(-y^) then the map 
has two fixed points x\ and x 2 on the interval (0, 1) with one on either side of the maximum, 
i.e., < X\ < exp(— < x 2 < 1. The map also has a fixed point at x = and while 
this is outside the domain (0,1) it is an attracting fixed point for x G (0,xi). Numerical 
underflow can be avoided by selecting parameter values 7 and (3 together with the initial 
point x(0) so that the twp conditions x\ < /(^) and x 1 < x(0) < ^ are satisfied. Future 
iterates will the be confined to x G (x±, ^). 

The function f(x) is displayed in Fig.l for 7 = 3.78 and several values of (3 = 0.2, 0.6, 1.0, 
1.5, 2.0. From this figure it can be seen that the main effect of increasing (3 at a fixed value 
of 7 is to shift the position of the maximum to higher values of x. In contrast to the logistic 
map, the Kim-Kong map is not symmetric about the maximum. In Fig. 2 we have displayed 
an empirical unimodal one-dimensional map that was obtained from Belousov-Zhabotinskii 
reaction experiments [12 — 14]. By comparing Figs. 1 and 2 it can be seen that the functional 
form of the Kim-Kong map for < (3 < 0.2 and 7 = 3.78 is similar to the experimentally 
obtained map. 

The Kim-Kong map undergoes a period doubling cascade to chaos at fixed (3 as 7 is 
increased. An example of the period doubling cascade is shown in in Fig. 2 of [6] for (3 = 0.2 
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and < 7 < 3.78. Denoting ^ k as the parameter at which the period 2 k 1 cycle becomes 
unstable we compute the sequence 

7fc+i - Ik 
Ok = 

lk+2 - lk+1 

from which we deduce that 5 = lim^oo 5k ~ 4.66920 in agreement with Feigenbaum's 
universal scaling result [5] for dissipative uni-modal maps. The other scaling exponent, a ~ 
2.50290 for successive branching widths is also recovered in the Kim-Kong map. 

We now consider a multifractal analysis[7 — 9] of chaotic orbits of the Kim-Kong map. 
First we review the definition of generalized dimensions in the multifractal formalism. Sup- 
pose that we divide the unit interval into M(e) cells of size e = j^- For a given chaotic orbit 
of length N let Pi = denote the fractional number of points in the ith cell, [(i — l)e, (i)e\. 
The generalized dimensions are now defined by 

D q = \im . g ^=] P \ (3) 
q e-o (q - 1) log e 

and the multifractal measures f q and a q are related to D q via the Legendre transform: 

f q = q^-l(q-l)D q ]-(q-l)D q (4) 

and 

<* q = -^[{q-l)D q ]. (5) 

We have made use of Eqs.(3)-(5) to compute the multifractal measures for the Kim-Kong 
map and the fractal dimension, D , which is numerically compared with that of other maps. 

In our numerical computations we have restricted ourselves to three cases for the two 
control parameters in the chaotic regime; (a) 7 = 3.0 and (3 = 0.12, (b) 7 = 3.4 and (5 = 0.12, 
and (c) 7 = 3.5 and f3 = 0.12. The multifractal measures that we have computed are based 
on 2 x 10 4 iterations of the map in Eq.(l). Figs. 4 and 5 are, respectively, the plots of 
D q versus q and the spectrum f q versus a q for chaotic orbits in the Kim-Kong map. The 
fractal dimension D or f ( i.e., the maximum value of D q or f q ) of the Kim- Kong map is 
numerically compared with that of other maps, in Table 1. 
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In conclusion, we have computed chaotic and multifractal properties of the Kim-Kong 
map - a two-parameter map of the unit interval onto itself. Although this map was investi- 
gated as a mathematical example in this work the map could prove useful as a model for the 
transport of quantum excitations, directed polymer problems, and the electron localization 
in quantum mechanics[18 — 20]. In future work on the Kim-Kong map we plan to investigate 
chaotic orbits of this map as examples of deterministic diffusion including the possibility 
of biased diffusion and anomalous diffusion. The chaotic orbits of the Kim-Kong maps 
will also be compared with randomness in a reaction-diffusion system with multi-species 
reactants[21 — 22]. 
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FIGURE CAPTIONS 



Fig. 1. A plot of x(n+l) verse x(n) graph in the Kim-Kong map [6], x(n+l) = 7exp[— f3 (log x(n)) 2 ][l — 
exp[-/3(logx(n)) 2 ]], for 7 = 3.78 and several values of (3 = 0.2, 0.6, 1.0, 1.5, 2.0. 

Fig. 2. Plots of two types of maps for chaotic states with impure malonic acid in the 
Belousov-Zhabotinskii reaction[12], where the circle and triangle data are the different U 
sequence of periodic states, respectively. 

Fig. 3. Plots of D q and q for 7 = 3.0 and f3 = 0.12 (thick solid line), 7 = 3.4 and (3 = 0.12 
(dashed line), and 7 = 3.5 and (3 — 0.12 (thick dashed line). 

Fig. 4. Plots of f q versus a q obtained from the chaotic sequences of the Kim-Kong map for 
three cases of the parameters; (a) 7 = 3.0 and f3 = 0.12 (thick solid line), (b) 7 = 3.4 and 
(3 — 0.12 (dashed line), and (c) 7 = 3.5 and /3 = 0.12 (thick dashed line). 



TABLE CAPTIONS 



Table 1. Summary of values of the fractal dimensions of chaotic orbits in the Kim-Kong map 
and other maps. 
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Maps Control Parameter Fractal Dimension 



Kim-Kong Map 


(a) 7 


= 3.0,13 = 0.12 


0.8395 




(b) 7 


= 3A,(3 = 0.12 


0.9790 




(c) 7 


= 3.5,(3 = 0.12 


0.9951 


Henon map [15] 


a = 


= 1.4,6 = 0.3 


1.26 


Kaplan- Yorke Map [15] 




a = 0.2 


1.431 


Logistic Map [16] 


b = 


3.5699456... 


0.538 


Lorentz Equation[15, 17] 






2.06 ±0.01 


Zaslavskii Map [16] 






1.39 
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